Fermi Liquid instabilities in two-dimensional lattice models 
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We develop a procedure for detecting Fermi liquid instabilities by extending the analysis of Pomer- 
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1. INTRODUCTION 

The Landau theory of the Fermi Liquid (FL) is one of 
the most important frameworks to understand conven- 
tional weakly interacting metallic systems^. The low en- 
ergy physics of interacting fermions in three dimensions is 
usually described by Landau's FL theory whose central 
assumption is the existence of single particle fermionic 
excitations, or "quasiparticlcs" , with a long lifetime at 
very low energies. In lower dimensions, however, the 
situation is much more interesting: in one dimensional 
systems Landau's quasiparticles are typically unstable, 
giving rise to the so called Luttinger Liquid. On the 
other hand, two dimensional lattice models are far more 
complicated to treat since conventional perturbation the- 
ory breaks down for densities close to half-filling, where 
competing infrared divergences appear as a consequence 
of Fermi Surface (FS) nesting and van Hove singularities. 

In Refi^ Pomeranchuk developed a method to diag- 
nose instabilities of the FL, by "deforming" the FS and 
studying the resulting energy gain. In its original form, it 
applies to systems with a three dimensional spherical FS, 
but it can be easily generalized to the two-dimensional 
continuum case. 

The experimental observation of exotic phases in 
strongly correlated systems has triggered an enormous 
effort from the theoretical community to try to under- 
stand their microscopical origin. One possible route to 
detect instabilities of a FL is precisely the analysis done 
by Pomeranchuk. Due to that, the Pomeranchuk insta- 
bility has been studied by several authors with different 
techniques in the last few years 3 ^^^^ and in particu- 
lar, the instability of the FL towards the nematic phase 
was investigated for several model s 10 ^ 1 ' 12 ! 13 ! 14 ! 15 . 

In this paper we develop a general method to trace 
such instabilities in lattice models, in a simple and rig- 
orous way. It allows for the study of systems which have 
an arbitrary shape of the FS in the absence of interac- 
tions, thus being applicable to models relevant to high 
temperature superconductors, manganites, ruthenates, 
etc . 16 d 7 i 18 i 19 . as long as one can rely on a perturba- 
tive analysis. It can be applied in principle to any lattice 



problem in a systematic way. 

We test our method within two examples, the attrac- 
tive Hubbard model and a model with forward scattering 
interactions that give rise to <i-wave FS deformation (the 
so-called "<i-wave Pomeranchuk instability"). 

The paper is organized as follows: Section [2] contains a 
detailed derivation of the method, with the proof of our 
formulas in subsection l2.11 and a shorthand recipe for the 
application of the results in subsection l2.2l Then in Sec- 
tion [3] wc apply the method to a two dimensional square 
lattice with the various interactions studied ir*i£, the s- 
wave interaction being studied in subsection 13. 2[ while 
the d-wave instability in subsection 13.31 Finally, section 
H] contain the conclusions, and some specific calculations 
are presented in the Appendix. 



2. TWO DIMENSIONAL POMERANCHUK 
INSTABILITY 

2.1. Derivation of the method 

In the theory of Landau's FL, the free dynamics at 
zero temperature is determined by the dispersion relation 
e(k). In terms of it, the FS is defined as the set of points 
in momentum space satisfying the equation 



H-e{k) =0. 



(2.1) 



In the ground state of the system, all single-particle states 
inside the FS /i — e(k) > are occupied, while those 
outside FS ^ — e(k) < are not. Excited states of the 
system are built by moving some particles from the inner 
single-particle states to the outer ones. 

The energy of such excited state as a functional of the 
change in the equilibrium distribution function can be 
written as 

E= d 2 k(e(k)-^)5n(k) + \- 



d\\d 2 k' f(k,k') Sn(k)Sn{k'), 
(2.2) 



where 8n(k) is the change in the distribution function 
n(k), and we have assumed that only two-particle inter- 
actions are present. The interaction function f(k, k') can 



be related to the low energy limit of the two particle ver- 
tex. 

Pomeranchuk criterion allows to identify low energy 
excited states of the system that make (|2.2p negative. 
This signals an instability, and the breakdown of the FL 
description. In what follows, we will carefully go through 
all the steps needed to perform such analysis. 

First let us define, associated to any given state of the 
system, a smooth function g(k) such that it takes posi- 
tive values at occupied single-particle states and negative 
values at unoccupied ones. Then at the frontier between 
these two regions we have the equation 



g(k) = 0. 



(2.3) 



For the ground state, such frontier coincides with the FS 
allowing us to choose 



g{k) =n- e(k) 



(2.4) 



Under a variation 5n{k) of the distribution function 
n(k), we get an excited state that can be described in 
terms of a new function g'(k) = g(k)+8g(k). The frontier 
between occupied and unoccupied single-particle states is 
now located at points satisfying 



g'(k)=g(k) + Sg(k) = 0. 



(2.5) 



By an abuse of language we will call the solution of this 
equation the deformed FS. 

Since at T — we have 5n{k) = ±1 we can write 



Sn(k) = H[g\k)] - H[g(k)] 



(2.6) 



where H (x) is the unit step function, defined by H(x) = 1 
if x > and H (x) = if x < 0. This can be replaced 
in (|2.2|1 to write the energy of the quasiparticles as a 
functional of g(k) and g'(k), namely 



E = J d 2 k (s(k) ») (H[g\k)} H[g(k)}) + ljd 2 kj d 2 k'f(k, k')(H[g'(k)] H[g(k)]){H[g\k')] H[g(k')\). (2.7) 



r 



To go further, we have to take into account the constraint 
imposed by the Luttinger theorem^, or in other words 
the preservation of the area of the FS under the defor- 
mation 



d 2 k5n{k) = 0. 



(2.8) 



By using (|2.6| this can be rewritten as a functional con- 
straint on the functions g{k) and g'(k) 



d 2 k H[g'(k)] = / d 2 k H[g(k)] 



(2.9) 



In two-dimensions the constraint (|2.9[) can be easily 
solved as follows. We first rename the integration vari- 
ables on the right hand side to k'. Next we assume that 
a change of variables k' = k + Sk(k) can take the right 
hand side into the form of the left hand side. Writing 
g'(k) = g(k) + Sg(k) we get two unknown functions to 
be solved for, namely Sg(k) and 5k(k), together with the 
equation 

d 2 kH[g(k)+Sg(k)] = d 2 k\l + djStf] H[g(k + Sk(k))] , 

(2.10) 

where i,j £ {1, 2} label the orthogonal directions in mo- 
mentum space. 

A particular class of solutions can then be obtained by 
solving the following equations 



ll + diSVl 



1 



g(k) + Sg(k) = g(k + Sk(k)). 



(2.11) 



The first line (|2.11|) implies that the change of variables 
going from k' to k is an area preserving diffcomorphism. 
The second line on the other hand, can be interpreted as 
saying that the variation Sg(k) is a translation of g(k) by 
an amount 5k. We can solve (J2.11J) as 



SW = (e e ' kd ' xdk - l)k l , 
Sg - {e^ d ^ -l)g. 



(2.12) 



Where A is a free function parameterizing the deforma- 
tion. If we assume that the deformation of the FS is 
small, then Sg(k) is also small and we can parameterize 
it in terms of a slowly varying A 



8k l ~ e v djX, 
Sg ~ e^djXdig. 



(2.13) 



In what follows, each specific form of A will characterize 
an excited state, the sign of the resulting energy will give 
us information about the instabilities. 

Now that we have solved the constraint, we go back 
to the energy of the quasiparticles (|2.7[) and write it in 
terms of the free unconstrained variable A. To simplify 
the resulting expression, we need to change variables to a 
more convenient coordinate system in momentum space. 
We choose a special set of variables 



g = g(k x ,k y ). 

■S = S(K X , Ky) . 



(2.14) 



where the new variable g varies in the direction trans- 
verse to the unperturbed FS. The variable s varies in 
the longitudinal direction tangent to the FS, namely it 
satisfies disdig = 0. 

Separating the energy (|2.7[) into a linear and an inter- 
action term E = L + I, we get for the linear part 

L = fd 2 k (e(k) - n) (H[g + Sg] - H[g\) = 



ds dg J(a, g){e{s, g) - y) (H[g + 5g] - H[g}) 



= \ ds \ dgj(a,g)(s(s,g)-fi) 

-Sg 



(2.15) 



where J = \d(k x , k y )/d(g, s)\ is the Jacobian of the trans- 
formation (|2.14p . Expanding in powers of the integration 
variable g around the unperturbed FS g = we get 

L= fdsf dgd- g [J( Sl g)(e(s,g)-f,)}- g= og + 0(6g 2 ) = 

J J-Sq 



ds[J( S ,g)Sg 2 } g = g ,= Q + 0(6g 3 ) 



(2.16) 



where in the second line we have integrated out the vari- 
able g and made use of the fact that (e(s, g = 0) — /i) =0. 
In order to replace the explicit form of 5g eq. (|2.13|) in the 
integrand of (|2.16|) we make use of the identity 

e lJ digdj\ = e lJ (digd g g + disd s g)(djgd g X + djsd s X) = 



e^digdjsd s \ = J 1 d s X. 



3 

(2.17) 



where we have used the fact that, according to our def- 
initions, d g g — 1 and d s g = 0. Now replacing in (|2.16D 

we get 



L 



ds[J- 1 (g,s)(d s X) 2 ] (2.18) 



The calculus of / is analogous and gives 



/= - J ds I ds' {f(g,s;g',s') (d s X)(d 3 ,X)]\ g=g , =0 .(2.W) 



Adding the two contributions we finally have 



E 



l - ds ds (f(0,s;0,s') + J- 1 (0,s)8(s-s')) x 



xd e X(0,s)d s ,X(0,s') 



(2.20) 



As the functions A(0, s) characterizing the excited states 
are arbitrary, we can equally work with functions ^>(s) = 
d s X(g,s)\ g =Q. In what follows we will be interested in 
excited states such that ip(s) £ L 2 [0, S]. Assuming that 
s makes a complete turn around the FS when it runs 
from to S, we also need to impose periodicity in that 
interval. 

Since the sign of E in cq. (|2.20|) determines the sta- 
bility of the FL, from all the above we conclude that the 
stability condition reads 



ds' j 'd S ^(s')^Lj- 1 (s)5{s - s') + /(«,«')) ^(.s) > 



(2.21) 



r 



where have we defined 



/(«,«') = f(9,s;g',s')\ g=gl=0 



1 9=0 



(2.22) 



Note that the stability condition has two terms, the 
first of which contains the information about the form 
of the FS via J _1 (s), while the second encodes the spe- 
cific form of the interaction in /(s, s'). There is a clear 
competition between the interaction function in the sec- 
ond term of the integrand and the first term that only 
depends of the geometry of the unperturbed FS. 

We see that E is a bilinear form, acting on the real 
functions tp(s) parameterizing the deformations of the 
FS 



E=(il>,il>), 



(2.23) 



where 



(u, v) — ds' dsu(s)- l/(s,s') + J \s)S(s — s')\ v(s') . 

(2.24) 



The stability condition is then equivalent to asking this 
form to be positive definite for any possible deformation, 
i.e. 



Vip :(i/),ip)>0. 



(2.25) 



In consequence, the natural way to diagnose an instabil- 
ity is to diagonalize this bilinear form and to look for 
negative eigenvalues. 

We can expand the functions tp(s) in some basis of 
i 2 [0, S] that we will denote {&(s)} 



i>(s) = Yl °i&( S ) 



and then write 



E= X] a *i a *2 (&l'&: 



(2.26) 



(2.27) 



The bilinear form ( , ) can be taken as a pseudo-scalar 
product, which is linear and symmetric but, in general, 



not positive-definite. Only in the free case f(s,s') = 
the positivity is ensured. If {£i(s)} are taken to be 
orthogonal with respect to this pseudo-scalar product, 
then the functional (|2.2ip is given by 



E 



E 



2 fJ, 
a i Xi 



(2.28) 



where xt = (Ch £«} is the square pseudo-norm of the or- 
thogonal functions. If Xi nas a negative value for some 
i, then by choosing the corresponding a 2 — 1 and a 2 = 
for j ^ i, we see that the energy is negative denoting an 
instability. In this case we say that we have an insta- 
bility in the i-th channel. In other words, the stability 
condition has been mapped into 



V*:xf >0, 



(2.29) 



the Xi being taken as the stability parameters. If any of 
these quantities is negative, then the FS is unstable. 

We perform such diagonalization by choosing a basis 
on Li [0, S] as a given set of functions {ipi} and then mak- 
ing use of the Gram-Schmidt orthogonalization proce- 
dure to transform it into an orthogonal basis {&}■ Note 
that, being the bilinear form not necessarily positive def- 
inite, the new basis cannot be normalized to 1 but to 
±1. 

This is our main result. Our method to search for 
Pomcranchuk instabilities, can be summarized in the fol- 
lowing recipe: 



2.2. Recipe 

1. Get the dispersion relation e(k) and the interaction 
function /(fc, k') for the model under study. 



2. Change variables according to ()2.14|) . The variable 
g is completely fixed by the dispersion relation ac- 
cording to (|2.4|) . The choice of s is arbitrary ex- 
cept for the constraint of being tangential to the 
FS, d iS d i9 = 0. 

3. Write the bilinear form E as in (|2.21|) . 

4. Choose an arbitrary basis of functions {ipi} of 
L 2 [0,5]. 

5. Apply the Gram-Schmidt orthogonalization proce- 
dure, verifying at each step whether (£j,£i) > 

6. If for a given channel i one finds that (£i,£j) < 0, 
the FS is diagnosed to be unstable. 

Note that since L2P, S] is infinite dimensional, the 
present method is not efficient to verify stability: at any 
step i it may always be the case that for some j, xt+j < 0- 
Moreover, we have not exhausted all the possible solu- 
tions of the constraint (|2.9p but only explored a subset 
of them. 



3. INSTABILITIES IN THE SQUARE LATTICE 

3.1. Contribution of the free Hamiltonian 

We start considering free fermions in the square lattice, 
with a Hamiltonian given by 



where 



H Q = YVe(fc) - v)c\c k . 



e(k) = — 2£(cos k x + cos k v ) , 



(3.1) 



(3.2) 



where only hopping to nearest neighbors has been taken 
into account. The FS is defined by 

g(k) = n — e{k) = fx + 2t(cosk x + cosk y ) = 0, (3.3) 

where fj, is the chemical potential. Notice that g > 
inside the area bounded by the FS, negative outside it, 
and zero at the FS. 

Now we follow the recipe given in Scction r2.2l changing 
variables according to (|2.14[) 

g(k x ,k y ) = fi + 2t(cos k x + cos k y ) , (3.4) 

/tan(fcj,/2)~ 



S\fV X: Ky) 



arctan 



\tan(fc x /2) 



(3.5) 



It is straightforward to see that g and s are mutually 
orthogonal variables. Using the following shorthand no- 
tation 



a 



cos k x , 
cos A: 



!) • 



(3.6) 
(3.7) 



we can write 



g = n + 2t{a + f3), 

and the Jacobian takes the form 

a/3-1 



J = t 



a 2 +/3 2 



(3.8) 



(3.9) 



Notice that J > 0. Writing a and (3 as functions g and s 
we have for the Jacobian evaluated at g = 0. 



J[g = 0,s} = 



1 



2t v /l-/3(/i) cos 2 (2s) 



(3.10) 



with /3(/x) = 1 — (^j) 2 - The limits for the variable s can 
be taken as — 7r < s < 7r. 

The inverse of the jacobian J -1 ^) can be expanded 
in series of sin(ns) and cos(ns) and it is straightforward 
to show that only the coefficients of cos(4ns) are non- 
vanishing. This results in the following expansion 



J- 1 (s)=^^cos(4n S ) ; 



(3.11) 
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FIG. 2: Phase diagram for f(s,s') = Uo, are displayed the 
regions of instability for the first three channels. For smaller 
interactions, higher channels are unstable closer to half filling. 



FIG. 1: The instability parameters. For Uo — we show the 
first 10 parameters xt as a function of /j,. For other values 
of the interaction we show only the parameters corresponding 
to the three lower channels that show instabilities, namely 
channels xo, X& an d Xi6 (colors are identified in Fig[2}. Notice 
that when we increase Co, the FL breakdown occurs first for 
the higher channels and closer to half filling. 



where the coefficients j% are fixed by the expansion. 
Some of them are given in the Appendix. The simplicity 
of this expansion suggest to take as our starting base in 
the Gram-Schmidt orthogonalization procedure the set 
{sin(ns),cos(ns)}. 

In the next subsections we will analyze as an exam- 
ple of application the possible instabilities in this two- 
dimensional fermion model when subjected to various in- 
teractions. In particular we are interested in interactions 
of the formic. 



f(k,k')= Constant x d{k)d(k') , 



(3.12) 



with 



d(k) = 1 , 

d(k) = (cos k x + cos k y ) . 

d(k) = (cosk x — cos k y ) , 



(s-wave) , (3.13) 

(extended s-wave) , 
(d-wave d x 2_ y 2) . 



3.2. s-wave instability 

First we consider a constant interaction corresponding 
to take^ d(k) = 1 so that 



f(k,k') = U , 



(3.14) 



Uq being a constant measuring the strength of the inter- 
action. This form of the interaction can be obtained by 



a Mean Field approximation or a first order perturbative 
expansion for the interaction function on the Hubbard 
mode l 21 i 22 . 

Using the product defined in eq. (|2.24[) we can calcu- 
late the first instability parameters. 



Xl,2 



u 

X3 



II 

Xi 



LL 

^5,6 



\7 



xt 



etc. 



27r(17o7r + #) 



(3.15) 










2(7rC/o+J M ) 2 vUo+ji; 



-iff 



./i2 
*"Ji 



TJ2 



2(wU +j£y 



(3.16) 



The stability parameters for the first channels are 
shown in Fig [TJ For Uq = all the Xn are positive as 
expected. When we increase the interaction, among the 
first 20 parameters, only Xqi X$ an( i Xi6 change. For 
simplicity only these parameters are plotted for Uq 7^ 0. 

With these first instability channels we can draw a 
qualitative phase diagram in the (/i, Uq) space as in Fig[2] 
where the first instability zones are shown and a tentative 
global phase diagram is drawn. 

Note that, when the interaction is increased, the first 
instability channel corresponds to the highest of the three 
shown in the figure {e.g. X\r)- This behavior is main- 
tained for channels xt with higher index z, and we as- 
sume that generically these higher channels will show the 



instability closer to half filling and for interactions arbi- 
trarily small. On the other hand, the higher the channel, 
the closer the instability region is to half-filling. Extrap- 
olating this behavior we see that the instabilities on the 
large- z channels take place only very close to \i = 0. 

The behavior for the extended s-wave with a form fac- 
tor d(k) = (cosk x +cosk y ) = a -\-/3 can be studied writing 
d(k) in terms of the variables g and s using the solutions 
of (|3.8p and evaluating at g = 0. We have d(s) = -^ and 
the interaction function reads 



/(..I -«■(£)' 



(3.17) 



Again f(s, s') is independent of the variables s and s' but 
now dependent of the chemical potential fi. 

The corresponding instability parameters are obtained 

by changing U — > U (^) in (|3.15p . and the phase dia- 
gram can then be inferred to be analogous to that of Fig. 

[2] but with the vertical axes replaced by Uq (£) . 



3.3. d-wave Pomeranchuk instability 

Now we investigate <i-wave Fermi Surface Deformation 
(dFSD) instability^ in the charge channel on a square 
lattice. The forward scattering interaction driving the 
dFSD has the formic 



f(k,k') = -Gd(k)d(k') 



(3.18) 



with G > and <i-wave form factors d{k) = (cosk x — 
cos k y ). The above expression for this effective interac- 
tion was obtained by Mctzncr et a&^ using functional 
renormalization group methods. 

Using the shorthand notation (|3.6p the interaction 
reads 



f(s,s') = -G(a-/3)(a'-/3') 
and using the solution of eq. (|3.8p we have 
2 fJ- l {s) 



d(s) = 



cos(2s) 



-1 



(3.19) 



(3.20) 



Notice that this interaction contains the Jacobian but its 
origin is totally independent of the treatment developed 
in the last sections. 

The form factor d(s) can be expanded in a series of the 
form 



d(s) 



OO 

E 

k=0 



d k cos((Ak + 2)s) , 



(3.21) 



where the first coefficients of the expansion are presented 
in the Appendix. 

Performing the Gram-Schmidt orthogonalization as in 
previous case, we find the x A '"P araine ters corresponding 
to this interaction. The results arc very similar, and we 




n/t 



FIG. 3: Phase diagram for the d-wave interaction. The first 
three unstable channels are shown. The dashed line corre- 
sponds to the the critical value for the interaction parameter 
studied in reference [ill within a Mean Field treatment. 



will only display here the first two parameters that show 
an instability of the system, namely 






Li 

Vs 



2GVd3+7TJ + 



TJl 



(3.22) 



(ji 2 -4jo 2 + 2jiJ2+ 



2gnM 2 -2(2j +j 1 ) 

+ J2 2 - jiJ3 - 2j (ji + :h) + ffTr (2di 2 (2j + h) 
- 8d d 1 {n +j 2 ) + W (2 j +j 3 ))) (3.23) 

Again, by making use of this first instabilities we can 
sketch the phase diagram corresponding to this interac- 
tion, as shown in Fig. [3] Similarly to the previous case, 
as the interaction grows instabilities appear first in the 
higher channels. The dashed line corresponds to the crit- 
ical value of the interaction found in ref [ll| by means of 
a Mean Field procedure. Notice that the lowest channels 
cover most of the instability zone. The phase diagram 
shown in Fig. [3] is consistent with the results presented 

in [mm. 

Unlike treatments using Mean Field, with the present 
formalism it is possible to identify the region in the pa- 
rameter space where each channel presents a breakdown 
of the Fermi liquid behavior. 



4. SUMMARY AND CONCLUSIONS 

In this paper we have developed a general procedure 
for detecting instabilities in two dimensional lattice mod- 
els. It is an extension of the formalism of Landau- 
Pomeranchuk, in particular for lattice systems with an 
arbitrary shaped FS and allows to describe the phase 
diagram of the system as an alternative to the usual pro- 
cedures. The steps are simple and applicable to a wide 
variety of systems. 

Complcmcntarily to other descriptions, our procedure 



permits to identify the breakdown of the Fermi liquid 
behavior on each instability channel independently. 

As a form of testing our procedure, we have analyzed 
the stability of the Fermi Liquid in a square lattice, for 
various interactions already studied in the literature. The 
s-wave instability in the electronic channel and the insta- 
bility produced by d-wave forward scattering interactions 
were studied at T = 0. The instabilities corresponding to 
low channels produce a breakdown of the FL behavior for 
a wide range of fillings, while those occurring for higher 
channels are closer to half filling. Our result are in good 
agreement with those obtained by different methods that 
were previously published by other authors. 

Generalization to higher dimensions, spin-dependent 
interactions or finite temperature can be achieved fol- 
lowing the same lines and the results will be presented 
elsewhere 2 ^. 



5. APPENDIX: ORTHOGONAL BASIS AND 
SERIES EXPANSION. 

In this Section we present the coefficients in the expan- 
sion of the functions used in this paper. For the Jacobian, 
the series takes the form 



J 1 ( S )=J2 j£cos(4ns), 



(5.1) 



where the first coefficients in the expansion are given by 



Jo 






1e 

7T 


i - 


n 2 ~ 

16 










±(me 


1- 


16" 


-4E 


['-£1 

16 


-H 


n i 


Fi 


1 
~2' 


5,2,1- 

2' ' 


16' 


+ 



2tt 2 Fi 



-1*2 1-*- 

2' 2' ' 16 






- | 

7T 



16 



47T 2 Fl 

3tt 2 Fi 



1 3 m 2 

9 1 — - 

2'2' ' 16 

_I 5 /^ 

2'2' ' 16 






— [E 

7T 



1-0- 

16 



+4tt [ 9 2 Fi 
- 18 a Fi 
+ 10 2 Fi 



■i.5,2,1-^ 

2'2' ' 16 



-i.5,.,i-^ 

2' 2' ' 16 

-11,4 1-^ 
2'2' ' 16 



(5.2) 



where E [m] is the complete elliptic integral 

E [m] = J yd- to sin 2 (f) dt , (5.3) 



and 2 Fi(a, 6; c; z) is the hypergeometric function 

Tic] 
a Fi(a,&;c;z) = 



(5.4) 



r(6)r(-6 + c) 

x /" i 6 - 1 (l-i) c - 6 - 1 (l-te)- a dt. 
./o 



The form factor for the d-wave forward scattering in- 
teraction can be expanded as follows 



d(s) = ^d k cos((4fc + 2)s) , 



(5.5) 



with 
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The orthogonal basis {^} depends of course on the spe- 
cific form of the interaction, but in all the cases studied 
here it satisfy the following properties: 

1 - The functions are either linear combinations of 
sin(s) or of cos(s) separately. There is not mixtures of 
sin and cos. 

2 - All the functions reduce to the expressions corre- 
sponding to the free case in the limit when the interaction 
parameter is sent to zero. 
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